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ABSTRACT 

We investigate the feasibility of detecting 21 cm absorption features in the afterglow spectra of 
high redshift long Gamma Ray Bursts (GRBs). This is done employing simulations of cosmic 
reionization, together with estimates of the GRB radio afterglow flux and the instrumental 
characteristics of the LOw Frequency ARray (LOFAR). We find that absorption features could 
be marginally (with a S/N larger than a few) detected by LOFAR at z L 7 if the GRB is a 
highly energetic event originating from Pop III stars, while the detection would be easier if 
the noise were reduced by one order of magnitude, i.e. similar to what is expected for the first 
phase of the Square Kilometer Array (SKAl-low). On the other hand, more standard GRBs 
are too dim to be detected even with ten times the sensitivity of SKAl-low, and only in the 
most optimistic case can a S/N larger than a few be reached atz~9. 
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1 INTRODUCTION 

Present and pla n ned radio facilities such as lofafQ 
( Ivan Haarlem et al.l l2013h . MW^I PAPElB and SKjfl aim 
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at detecting the 21 cm signal from the Epoch of Reioniza- 
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), or absorption features in the spectra of high-z radio- 


loud sources (e.g. Carilli et alj 

2 OO 2 I: Furlanettol 20061: Xu et al.l 

2 OO 9 I: Mack&Wvithe 

2 OI 2 I: 

MeiksinI 

20 111: Xuetal 201 ll; 

Vasiliev & ShchekinovI 

2012l;ICiardi et alj 

2013[;lEwall-Wiceetal.l 

2014[). These observations would offer unique informafion on the 


statistical properties of the EoR (such as e.g. its duration), the 
amount of Hi present in the intergalactic medium (IGM), and, 
ultimately, the history of cosmic reionization and the properties of 
its sources. 

In particular, the detection of absorption features in 21 cm 
could be used to gain information on the cold, neutral hydrogen 
present along the line of sight (LOS) towards e.g. high-z quasars, 
similarly to what is pres ently done at lo wer redshift with the Lya 
forest (for a review see lMeiksinll2009l) . While the detection and 
analysis of the 21 cm forest is in principle an easier task compared 
to imaging or even statistical measurements and could probe larger 
A:-modes, these absorption like experiments are rendered less likely 
to happen due to the apparent lack of high-z radio-loud sources (e.g. 
ICarilli et al.l l2002|jXti^talj2()09|). the most distant being TN0924- 
2201 at z=5d9 l lvan Breugel et al.ll 19991) . 

Gamma Ray Bursts (GR Bs) have been suggested as alterna¬ 
tive b a ckground sources ( e.g. Iloka & MeszarosIboOSi : llnoue et ^ 
l2007l: iToma et all |201l|). as they have been observed up 
to z ~ 8-9 S^aterra et al.l l2009l : iTanvir et alj 12001^ : 


ICucchiara et al.l 12011 ). t hey are expected to occur up to the 
epoch of the first stars tBrornrn_&lxteb| [2i)0j; iNatarajaneLalJ 


200j;_lKomissarov&^aAa^_[2010a|; |Meszaros_&Rees 
Camnisi et al .1 1201 ll: ISuwa & Iokall201 ll: iToma et alJl201 


ID 


20 id; 

and to 


be visible in the IR and radio up to very high-z due to cosmic time 
dilation effects jCiardi & LoeblH()od : lLamb & Reichartll200d) . Be¬ 
cause of the latter, if a GRB afterglow were observed in the IR by 
e.g. ALM^lat a post-burst observer time of a few hours, this would 
offer a few days to pla n for an observation in the radio accordingly 
(e.g. llnoue et alj2007l) . 

In this letter we investigate the feasibility of detecting the 
21 cm forest with LOFAR and SKA in the radio afterglow of high-z 
GRBs. In Section 2 we present the method used to calculate the for¬ 
est; in Section 3 we discuss the properties of the GRBs; in Section 4 
we present our results; and in Section 5 we give our conclusions. 


2 THE 21 CM FOREST 

In this work we make use of the pipeline developed in Ciardi et 
al. (2013, hereafter C2013) to assess the feasibility of detecting the 
21 cm forest in the spectra of high-z GRBs. For this reason, here 
we only give a brief overview of this pipeline, while we refer the 
reader to the original paper for more details. 

If a radio-loud source is located at redshift Zj, the photons 
emitted at a frequency v > V 2 icm = 1-42 GHz can be absorbed 
by the neutral hydrogen encountered along their LOS at z = 


^ http://www.almaobsei'vatory.org/ 


(V2icm/v)(l -I- z.i) - 1, with a probability (1 - The optical 

depth in the low-freq uency limit, as appropriate for the 21 cm li ne, 
can be written as (e.g. lMadau et all 1 9971 ; iFurlanetto et al.ll2006h : 


Tlic, 


3 h„dA2la 




32;r 


+ zXdm/druY 


( 1 ) 


where is the hydrogen number density, xhi is the mean neutral 
hydrogen fraction, Tj is the gas spin temperature, A 2 icm = 2.85 x 
10“'^ s“' is the Einstein coefficient of the transition, and dv\\ldr\\ 
is the gradient of the velocity along the LOS, with ry comoving 
distance and un proper velocity including the Hubble flow and the 
gas peculiar velocity. The other symbols have the usual meaning. 

The optical depth has been calculated using the simulation 
of reionization called 2)4.39 in C2013. This has been obtained by 
post-processing a GADGET-3 (an updated ve rsion of the publicly 
available code GADGET-2; see lSt)ringell2005h hydrodynamic sim¬ 
ulation with the 3D Monte Carlo radiative transfer code CRASH 


dCiardi et alJl200ll;lMaselli et al.ll200l l2009l; IPierleoni et aI1l2009l; 

IPartl et alJl201ll ; lGraziani et alj|201^ which follows the propaga¬ 

tion of UV photons and evaluates self-consistently the evolution 
of H I, He I, He ii and gas temperature. The hydrodynamic simu¬ 
lations were run in a box of size 4.39l!“* Mpc comoving, with 
2 X 256^ gas and dark matter particles, and cosmological param¬ 
eters Da = 0.74, D,„ = 0.26, D* = 0.024h^, h = 0.72, Uj = 0.95 
and <Tg = 0.85, where the symbols have the usual meaning. The gas 
density, temperature, peculiar velocity and halo masses, were grid- 
ded onto a uniform 128^ grid to be processed with CRASH. The 
sources are assumed to have a power-law spectrum with index 3. 
For more details on the choice of the paramete rs and the results of 
the reionization histories we refer the reader to lciardi et al.l ilQl'j) 
and C2013. Here we further note that the simulations were designed 
to mat ch WMAP observatio ns of the Thomson scattering optical 
depth dKomatsu et al.l 201 ih. while mor e recent Planck measure¬ 
ments dPlanck Collaboration et alJl2Q15h favor a lower value, i.e. 
a delayed reionization process. This does not change our conclu¬ 
sions, and it actually makes them conservative, as more H i would 
be expected at each redshift compared to the model considered 
here. 

Random LOS are cast through the simulation boxes and the 
corresponding 21 cm absorption is calculated. It should be noted 
that here we consider a case in which the temperature of the IGM 
is determined only by the effect of UV photons, i.e. gas which is 
not reached by ionizing photons remains cold, while the effect of 
Lyff and x-ray photons on the 21 cm forest is extensively discussed 
in C2013. 

Once the theoretical spectra are evaluated, instrumental effects 
and noise need to be included to calculate mock spectra. To do this 
we assume a radio source at redshift Zs, with a power-law spectrum 
with spectral index a, and a flux density S m(Zs)- We then simulate 
spectra as they would be seen by LOFAR. While we refer the reader 
to C2013 for more details, here we just mention that the noise <t„ 
is given b3B 




SEFD 


ns Y2N(N-l)uAv 


( 2 ) 


® We note that eq. |2]is con'ect as long as the S EFD is calculated theoreti- 
cally for a single polarization using the system temperature (as done in this 
paper), while when the S EFD is determined observationally from Stokes I 
the noise is reduced by a factor of sqrt(2), since it combines the two cross¬ 
dipole sensitivities. 
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where W ~ 1.3 is a factor which incorporates the effect of weight¬ 
ing, tis = 0.5 is the system efficiency, Av is the bandwidth, f,„, is the 
integration time, = 48 is the number of stations, and the system 
equivalent flux density is given by: 


SEFD = 




(3) 


where Kg is Boltzmann’s constant, Aeff = 1.5626) m? is 

the effective area of each dipol^il A^dip is the number of dipoles 
per station (we assume 24 tiles per station with 16 dipoles each), 
77a = 1 is the dipole efficiency, and the system noise is T^ys = [140 + 
60(v/300 MHz)“^'^^] k[|. As a reference, cr„ = 0.66 mjy for v = 
130 MHz, = 1000 h and Av = 10 kHz. 


3 RADIO AFTERGLOWS OF HIGH-Z GRBS 

Here we discuss aspects of t he radio afterglow emission from GRBs 
(see e.g. Ivan EeitM l2013l : iGranot & van der HorstI I20l4 for re¬ 
cent reviews) that are most relevant for studies of the 21 cm for¬ 
est at high redshift. GRB afterglows consist primarily of broad¬ 
band synchrotron emission from nonthermal electrons accelerated 
in the forward shock of relativistic blastwaves, which are driven 
into the ambient medium by transient, collimated jets triggered by 
the GRB central engine. The low-frequency radio flux is typically 
suppressed in the beginning due to synchrotron self-absorption and 
rises gradually as the blastwave expands. The light curve at a given 
frequency v peaks when the emission becomes optically thin to 
self-absorption, after which it decays, according to the overall be¬ 
havior of the decelerating blastwave. As background sources for 
observing the 21 cm forest, the emission near this peak flux time is 
naturally the most interesting. 

The expected observer peak time, fpt(v), and the correspond¬ 
ing peak flux, 5in(v), can b e evaluated u sing the formulation 
outlined in the Appendix of llnouel i2004h. which is ba s ed on 
standard discussions in th e literature (e.g. Isarfrt 'ZI |1998L[I^: 
IPanaitescu & Kurn^bOOOl) and is sufficient for our purposes of es¬ 
timating the radio afterglow emission at relatively late times after 
the burst. We consider cases for which the peak time occurs after 
the crossing time of the minimum frequency, v,„, as well as the jet 
break time but before the non-relativistic transition time, valid for 
the range of parameters chosen her^. For concreteness, the spectral 
index of the accelerated electron distribution is taken to be p = 2.2, 
implying a radio spectral index of a = 0.6 for the optically thin, 
slow cooling regime. We also choose =0.1 and eg = 0.01 for 
the fractions of post-shock energy imparted to accelerated electrons 
and magnetic fields, respectively, consiste nt with observationally 
inferred values l lPanaitescu & Kumaill200lh . This gives: 

t,,(v) = 540 days f (4) 



and 

5,n(v) ^ 1.1 X 10“^ mJy (5) 


^ Note that the same expression in C2013 contains a typo. 

^ We note that the values obtained from eq.[2]are very similar to the real 
ones reported in the LOFAR official webpage. 

^ Although the time evolution after the jet break i s not explicitly d escribed 
in llnoud i2QQ4l) . it is taken into account following ISari et aP il999h . 


Di^iz) 


3.3 X 1029 cm 


-2 / 1 , \ 2.53 

1 


11 


(lOS Hz) 


where E = 10^^ £53 erg is the isotropic-equivalent blastwave kinetic 
energy, 0 = 0.1 0 _i rad is the jet collimation half-angle, n,„egi,„„ = 
72o cm“^ is the ambient medium number density, and Di(z) is the 
luminosity distance. 

The known population of GRBs (referred to as GRBll, as they 
are expected to originate from standard Pop 11/1 stars) has been ob¬ 
servationally inferred to possess values of these quantities up to 
E ~ 10^** erg, 9 ~ 0.3 rad and/or down to ~ lO^’^cm^^ 

jPanaitescu & Kumail 1200 ll : iGhirlanda et al.l l2()13al) . Thus, they 
may provide fluxes up to 5in ~ 0.1 mJy at v ~100 MHz and 
tph ~ 3000 days (i.e. as sources virtually steady over several 
years) for events at z ~ 10 under favor able conditions (see also 
lloka & Meszarosll2005l : llnoue et ^|2007|) . 

On the other hand, although not yet confirmed by observa¬ 
tions, an intriguing possibility for high-z detections is the exis¬ 
tence of GRBs arising from Population 111 stars (GRBlll), with 
significantly larger blastwave energies, up to values as high as 
E ~ 10^^ erg, by virtue of their prolonged energy release fu¬ 
eled by accretion of the extensive envelopes of their progeni- 
tor stars 1 ^szaros & Reed 120101 : iKomissarov & BarkovI 1201 Obi : 
ISuwa & lokal 1201 ih . Compared to known GRBs, their blast- 
waves can expand to much larger radii and consequently with 
much brighter low-frequency radio emission, potentially exceed¬ 
ing S in ~ 10 mJy at v ~100 MHz and ~ 3 x 10"* days for events 
at z ~ 20 dToma et al.ll 201 ll : IGhirlanda et alj2013bh . 

Note that although more recent and detailed theoretical stud¬ 
ies of afterglow emission relying on hydrodynamical simulations 
have revealed n on-trivial deviations from the simple de scription 
presented above dvan Eertenll2013l : IGhirlanda et al.ll2013bl) . the lat¬ 
ter should still be sufficient for our aim of discussing expectations 
for observations of 21 cm forest. _ 

Finally, as reference numbers, ICampisi et alJ ( 1201 ih find that 
~ 1 (< 0.06) yr ' sr ' GRBll (GRBlll) are predicted to lie at z > 6 . 
This translates into ~ 7.5 x 10“^ GRBll (~ 4.5 x 10^^ GRBlll) per 
year in a 25 deg^ (LOFAR) field of view, and ~ 3 (~ 0.2) per year 
in a lO'' deg^ (SKA) field of view. 


4 RESULTS 

For an easier comparison to a case in which the background source 
is a radio-loud QSO, here we show results for the same LOS and at 
the same redshift of C2013 (see their Figs. 11, 12 and 13). 

In the upper panels of Figure[T]we plot the spectrum of a GR¬ 
Blll at Zs = 10 (i.e. v ~ 129 MHz) with a flux density 5in(Zs) = 
30 mJy. The simulated absorption spectrum. Sabs, is shown together 
with the observed spectrum. Sobs, calculated assuming an observa¬ 
tion time till, = 1000 h, a bandwidth Av = 10 kHz and a noise cr„ 
(given in eq. |2l le ft panels) and 0.1 cr„ (similar to the value expected 
for SKA1-1OW0. which is '-l/Sth of the LOFAR noise; right pan¬ 
els). In the lower panels of the Figure we also show the quantity 
l^in - Sabsl / l^obs “ ^absl. which effectively represents the signal- 
to-noise with which we would be able to detect the absorption. If 
indeed such powerful GRBs exist, then absorption features could 


SKAl-low is the first phase of the SKA covering the lowest frequency 
band. 
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129.6 129.8 130129.6 129.8 130 

[MHz] 


Figure 1. Upper panels: Spectrum of a GRBIII positioned at = 10 (i.e. 
V ~ 129 MHz), with a flux density i'inCzj) = 30 mJy. The red dotted lines 
refer to the intrinsic spectrum of the source, 5^', the blue dashed lines to 
the simulated spectrum for 21 cm absorption, ^abs^ the black solid lines 
to the spectrum for 21 cm absorption as it would be seen with a bandwidth 
Ay = 10 kHz, after an integration time tjni = 1000 h. The left and right 
panels refer to a case with the noise given in eq.[2|(LOFAR telescope) 
and with 0.\n (expected for SKAl-low), respectively. Lower panels: S/N 
corresponding to the upper panels. See text for further details. 



165.2 165.4 165.6 165.2 165.4 165.6 

^„b. [MHz] [MHz] 


Figure 2. As Figure [T] but the GRBIII is positioned at z, = 7.6 (i.e. v ~ 
165 MHz) and Av = 5 kHz. 

be detected by LOFAR with an averagl]]] signal-to-noise S/N~5, 
while if the noise were reduced by a factor of 10 the detection 
would be much easier. 

Figure shows a LOS to the same GRB located at z = 7.6, 
when the IGM is ~ 80% ionized by volume. The LOS has been 

We note that the definition of ‘average’ is somewhat arbitrary and de¬ 
pends on the frequency range used. Here the average refers to the one cal¬ 
culated over the frequency range shown in the Figures. 


GRBIII at z^=14 



^obs [MHz] [MHz] 

Figure 3. As Figure [T] but the GRBIII is positioned at Zs = 14 (i.e. v ~ 
95 MHz) and Av = 20 kHz. Note that the S/N in the lower right panel is 
always higher than the range covered by the axis. 

chosen to intercept a pocket of gas with Tzicm > 0.1 to show that 
strong absorption features could be detected, albeit by LOFAR only 
marginally with a S/N of a few, also at a redshift when most of the 
IGM is in a highly ionized state. 

On the other hand, a strong average absorption (rather than 
the absorption features seen in the previous figures) can be eas¬ 
ily detected (S/N > 10) in the spectra of GRBs located at very 
high redshift, as shown in Figure]^ If the intrinsic spectrum of the 
source could be inferred through other means, for example, accu¬ 
rate spectral measurements of the unabsorbed continuum at GHz 
frequencies and above, such detection could be used to infer the 
global amount of neutral hydrogen present in the IGM. 

We have applied the same pipeline also to more standard GR- 
BII, which have a flux density two to three orders of magnitude 
lower than a GRBIII. In this case we find that, even if we could 
collect 1000 h of observations with a noise 0.01 cr„ (i.e. 1/lOth of 
the SKAl-low noise), these would be barely enough to detect ab¬ 
sorption features in 21 cm. Also in the most optimistic case, with 
S in(z.i) = 0.1 mJy, a positive detection would be extremely difficult 
at any redshift, as shown in Figure|4| and only at z i 9 a S/N larger 
than a few could be reached. 


5 CONCLUSIONS 

In this letter we have discussed the feasibility of detecting 21 cm 
absorption features in the spectra of high redshift GRBs with LO¬ 
FAR and SKA. The distribution of H i, gas temperature and ve¬ 
locity field used to calculate the optical depth to the 21 cm line 
h ave been obtained f rom the simulations of reionization presented 
in iciardi et alj ( 1201 2h and C2013. We find that: 

• absorption features in the spectra of highly energetic GRBs 
from Pop III stars could be marginally (with a S/N larger than a 
few, depending on redshift) detected by LOFAR at z ~ 7; 

• the same features could be easily detected if the noise were 
reduced by one order of magnitude (similar to what is expected for 
SKAl-low); 


© 2014 RAS, MNRAS 000,[T]]6] 















































































































The 21 cm forest with LOFAR and SKA 5 



165.2 165.4 165.6 129.6 129.8 130 95.2 95.6 96 

[MHz] t;,,, [MHz] [MHz] 


Figure 4. Upper panels: Spectrum of a GRBII with a flux density S in(Zs) = 
0.1 mJy. The red dotted lines refer to the intrinsic spectrum of the source, 
5ini the blue dashed lines to the simulated spectrum for 21 cm absorption, 
^absi the black solid lines to the spectrum for 21 cm absorption as it 
would be seen after an integration time = 1000 h with a noise 0.01 (T„ 
(i.e. 1/lOth of the SKAl-low noise). The panels refer to a case with = 7.6 
and Ay = 5 kHz (left), = 10 and Av = 10 kHz (middle), = 14 and 
Av = 20 kHz (right). Lower panels: S/N corresponding to the upper panels. 
See text for further details. 


• the flux density of a more standard GRB is too low for ab¬ 
sorption features to be detected even with ten times the sensitivity 
of SKAl-low. Only in the most optimistic case with a flux density 
of 0.1 mJy, can a S/N larger than a few be reached at z 9. 


The problem of a small flux could be alleviated in case of a lensed 
GRB. Lensing of high-z sources has b een already discuss ed in the 
literature both from a theoretical fe.g. lWvithe et al.ll201lh. and an 
observational (e.g. with the Frontier Fields as in Qesch et alj20l4) 
perspective. Alternatively, a statistical detection of the 21 cm forest 
could be attempted, as already suggested by several authors (e.g. 

iMeiksinllioTIljMack & Wvithell201 IjEwall-Wice et alj2014h . 
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